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Abstract 

Background: Diapause is programmed developmental arrest coupled with the depression of metabolic activity and 
the enhancement of stress resistance. Pupal diapause is induced by environmental signals and is prepared during 
the prediapause phase. In the cotton bollworm, Helicoverpa armigera, the prediapause phase, which contains two 
sub-phases, diapause induction and preparation, occurs in the larval stage. Here, we performed parallel proteomic 
and metabolomic analyses on H. armigera larval hemolymph during the prediapause phase. 

Results: By two-dimensional electrophoresis, 37 proteins were shown to be differentially expressed in diapause-destined 
larvae. Of these proteins, 28 were successfully identified by MALDI-TOF/TOF mass spectrometry. Moreover, a total of 
22 altered metabolites were found in diapause-destined larval hemolymph by GC-MS analysis, and the levels of 17 
metabolites were elevated and 5 were decreased. 

Conclusions: The proteins and metabolites with significantly altered levels play different roles in diapause-destined 
larvae, including diapause induction, metabolic storage, immune response, stress tolerance, and others. Because 
hemolymph circulates through the whole body of an insect, these differences found in diapause-destined larvae most 
likely correspond to upstream endocrine signals and would further influence other organ/tissue activities to determine 
the insect's fact: diapause or development. 



Background 

The process of developmental arrest in insects, known 
as diapause, is a critical mechanism that allows for 
individual survival in harsh environmental conditions. 
Like hibernating mammals, diapausing insects can resist 
environmental stress by depressing metabolic activity and 
enhancing resistance to adversity. Diapause has evolved as 
a basic strategy for survival. Temperature, photoperiod, 
and nutrition are common environmental signals that 
induce diapause [1], Diapause is usually divided into 
three successive phases: prediapause, diapause, and 
postdiapause. The pre-diapause phase is further divided into 
two sub-phases: the induction phase and the preparation 
phase [2]. Environmental signals are perceived by insects in 
the prediapause phase, and physical signals are than 
transduced into humoral factors that trigger the decision to 
enter diapause. Individuals destined for diapause undergo 
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subsequent physiological and metabolic changes, including 
the storage of additional energy reserves, the deposition of 
extra layers for waterproofing the cuticle, and the synthesis 
of storage proteins in the hemolymph. Changes that take 
place during the prediapause phase are vital to individual 
survival. However, litde is known about the molecular 
mechanisms behind these physiological and metabolic 
shifts in the prediapause phase [3]. 

The cotton bollworm, Helicoverpa armigera (Noctuiadae, 
Lepidoptera), an agriculturally important pest, enters 
diapause at the pupal stage. Diapause is induced in 
this species by low temperature and short daylength 
during larval development. The prediapause phase in 
H. armigera occurs during the fifth and sixth instars: 
the induction phase bridges the fifth instar to the 
early sixth instar, and the preparation phase takes 
place during the mid-to-late stage of the sixth instar [4]. 
We have previously performed proteomic and metabolomic 
analyses of nondiapause- and diapause-destined larval 
brains from H. armigera. We found that proteins and 
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metabolites with significantly altered levels in diapause- 
destined individuals during the induction phase are most 
likely involved in a special memory mechanism used for 
the storage of photoperiodic information, and those 
during the preparation phase are most likely involved 
in the regulation of specific accumulation for energy 
reserves [5]. However, it is not yet known what occurs in 
the hemolymph as a result of these memory signals. 
In the present study, we focus on the changes in the 
proteomic and metabolomic profiles of the hemolymph in 
diapause-destined larvae of H. armigera in the prediapause 
phase. Hemolymph is present in the open circulatory 
system in insects to provide intermediary metabolites for 
individual growth and development, as all tissues or organs 
are exposed to hemolymph. The hemolymph also plays 
important roles in hormone and nutrient transport, energy 
deposition, innate immune responses, waste removal, and 
water balance. 

There are no obvious phenotypic differences to be 
observed between nondiapause- and diapause-destined 
larvae because both types of larvae have the same 
phenotypes of growth, development, and metamorphosis. 
Thus, it is difficult to identify genes that are specifically 
expressed or metabolites that are uniquely produced 
in diapause-destined larval hemolymph. Combined 
large-scale studies, such as transcriptomics, proteomics, 
and metabolomics, are particularly useful tools for the 
study of interactions between organisms and their 
environments. Recently, both proteomic and metabolomic 
techniques have been used successfully in the study of 
insect diapause [6-10]. Currently, however, little is known 
about the proteomic and metabolomic differences 
between hemolymph in nondiapause- and diapause- 
destined larvae in the prediapause phase. To understand 
the molecular events that occur in the prediapause phase, 
larval hemolymph from H. armigera was investigated by 
parallel proteomic and metabolomic analyses. These 
experiments were conducted on hemolymph from four 
stages of the prediapause phase: the mid-late stage of the 
fifth instar, the early stage of the sixth instar (belonging to 
the diapause induction phase), and the middle and the late 
stages of the sixth instar (belonging to the diapause 
preparation phase). A total of 37 differentially expressed 
proteins and 22 altered metabolites were identified, and 
these data provided valuable insights into the molecular 
characteristics of diapause induction and preparation. 

Results 

Two-DE analysis of differentially expressed proteins 

Two-dimensional electrophoresis (18 cm IPG strips, 
pH 4-7, linear) was combined with silver staining for in- 
vestigations of the proteomic differences in hemolymph 
at four larval stages: a) the mid-late stage of the fifth 
instar, b) the early stage of the sixth instar (diapause 



induction phase), c) the middle stage, and d) the late 
stage of the sixth instar (diapause preparation phase). 
Quantitative analysis revealed that 37 differentially 
expressed proteins were found in all three biological 
replicates in diapause-destined larval hemolymph during 
the prediapause phase; these proteins are marked with 
small arrows in the gel images (Figure 1). In the diapause 
induction phase, 2 protein spots were more prominent 
(spots no. al, a5) and 5 were less intense (spots no. a2-a4, 
a6, a7) in the mid-late stage of the fifth instar. Four 
protein spots were more intense (spots no. b4, b6-b8) 
in the early sixth instar, while 4 were less intense 
(spots no. bl-b3, b5). In the diapause preparation phase, 
7 protein spots were more pronounced (spots no. cl, 
c5-cl0) and 3 were less prominent (spots no. c2-c4) 
in the middle stage of the sixth instar. Seven protein 
spots were more intense (spots no. d3-d8, dll) in the 
late stage of the sixth instar, and 5 were less intense 
(spots no. dl, d2, d9, dlO, and dl2). Thus, a total of 
15 differentially expressed proteins were detected in 
the induction phase and 22 were detected in the 
preparation phase (Figure 2A). 

Protein identification and classification 

The 37 differentially expressed proteins were excised from 
the gels and identified through MALDI-TOF/TOF mass 
spectrometry. Representative PMF and MS/MS spectra of 
protein spot b3 are shown in Additional file 1: Figure SI. 
The MS/MS data were searched against an in-house 
H. armigera EST database constructed in our laboratory, 
the NCBInr database or the NCBI ESTothers database. 
A total of 17 proteins were successfully identified in the 
EST_Har database, and 11 proteins were identified in 
either the NCBInr database or the NCBI EST_others 
database (Table 1; Additional file 2: Figure S2). As shown 
in Table 1, several different spots in a gel were identified as 
the same protein. For example, both protein spots a3 and 
a4 were identified as glucosamine (N-acetyl)-6-sulfatase 
precursor, which is most likely a result of posttranslational 
modification. Using this method, a total of 28 differentially 
expressed proteins were identified in larval hemolymph. 
Because there are no commercial antibodies against these 
differentially expressed proteins from H. armigera, validation 
experiments such as western blot analysis were not 
performed in this study. 

The 28 identified proteins were classified according to 
the KEGG pathway maps or subordinately according to 
their biological processes from Gene Ontology. Proteins 
that were assigned functions, including glycan biosynthesis 
and metabolism, metabolism of cofactors and vitamins 
according to the KEGG pathway maps, or those related to 
metabolic processes according to Gene Ontology, were 
merged into a class referred to as "other metabolism". 
Among the differentially expressed proteins, 36% (10/28) 
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(See figure on previous page.) 

Figure 1 Representative 2-DE images of the proteins from H. armigera larval hemolymph. Proteins were extracted from nondiapause- (N) 
and diapause-destined (D) larval hemolymph at four stages: (A and B) mid-late stage of the fifth instar, (C and D) early stage of the sixth instar, 
(E and F) middle stage of the sixth instar, and (G and H) late stage of the sixth instar. Proteins were separated by IEF (18 cm IPG strips, pH 4-7, L) 
and 1 2% SDS-PAGE. Gels were stained by MS-compatible silver stain. Based on triplicate replications, only statistically significant (p-value < 0.05) 
protein spots that changed > 1 .5-fold were considered for further analysis. Differentially expressed protein spots are marked with letters and 
numbers. Letters 'a', 'b', 'c', and 'd' represent the mid-late stage of the fifth instar, the early, the middle, and the late stages of the sixth 
instar, respectively. 



were categorized as being involved in protein metabolic 
process. An additional 25% (7/28) were classified as being 
involved in metabolism of nutrients including carbohydrate 
metabolism, amino acid metabolism and other metabolism. 
Another 14% (4/28) were categorized as signaling molecules 
and interaction (Figure 2B). 

GC-MS analysis of significantly changed metabolites 

Comparative metabolomic analysis was also carried out 
between hemolymph samples from nondiapause- and 
diapause-destined H. armigera larvae in the prediapause 
phase. Metabolites were extracted from larval hemolymph 
at four stages, as described above, and analyzed by 
GC-MS. Representative GC-MS total ion chromatograms 



depicting the metabolite profiles of larval hemolymph are 
shown in Additional file 3: Figure S3. A total of 59 metabo- 
lites were identified across all four stages, including 18 
amino acids, 12 sugars and polyols, 6 metabolic intermedi- 
ates, 3 small molecules, and 20 other metabolites (Table 2). 
Based on four biological replicates, 22 metabolites exhibited 
significant changes in abundance in diapause-destined 
larvae. In the mid-late stage of the fifth instar, the 
only significant difference was found for phosphorylethano- 
lamine, as its content was increased in diapause-destined 
larvae. In the early stage of the sixth instar, 4 metabolites 
were detected at altered concentrations. Levels of trehalose, 
glucose, and mevalonic acid were higher, while less 
histidine was present. Thirteen metabolites were present 
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Figure 2 Number of differentially expressed proteins (A) and classification of the identified proteins (B) in H. armigera larval 
hemolymph in the prediapause phase. 
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Table 1 Differentially expressed proteins identified by MALDI-TOF/TOF 



Spot no. 


Protein name 


Organism 


NCBI 
accession no. 


Theoretical 
MW(kDa)/pl 


Fold 
change 


Database 


Mowse 
score 3 


No. of 
peptides 


EC no. 


Proteins increased in diapause larval hemolymph 


al 


Serpin 3a 


Manduca sexta 


gi|27733415 


51.13/526 


2.8 


ESTJHar 


311 


4 


- 


a5 


Ommochrome binding 
protein 


M. sexta 


—1 1 A r\r\/ — 7 n 

gi 400673 


30.61/5.89 


1 .5 


EST_Har 


369 


4 




b6 


Glutathione S-transferase 


Helicoverpa amigera 


gi]297340786 


23.72/6.18 


2.8 


EST_Har 


202 


3 


EC:2.5.l .1 8 


b7 


Hydroxypyruvate isomerase 


Danaus plexippus 


™ : 1 ~) / — 7 r f\ a a~7 r 

gi 357604475 


29. 1 2/6.00 


specific 


EST_others 


1 5 1 


2 


EC: 5. 3. 1 .22 


b8 


Pyridoxine 5'-phosphate 
oxidase 


Bombyx mori 


gi|l 1 2984194 




1 .5 


EST_Har 


120 


2 


EC:1 .4.3.5 


c5 


C-type lectin 7 


H. armigera 


gi|385202653 


34.99/5.98 


2.3 


ESTJHar 


116 


2 




c7 


Serpin 


H. armigera 


gi|314912147 


43.53/5.04 


2.7 


ESTJHar 


92 


2 




c9 


Imaginal disc growth 
factor-like protein 


Mamestra brassicae 


gi|85726208 


47.83/7.08 


3.8 


NCBI 


99 


1 b 


EC3.2.1.14 


d4 


Beta-1,3-glucan 
recognition protein 3 


H. armigera 


gi|208972535 


53.37/5.85 


1.9 


NCBI 


262 


4 


EC3.2.1- 


d5 


Putative serine protease-like 
protein 2 


Lonomia obliqua 


gi|56462300 


32.15/5.79 


5.9 


EST_others 


182 


2 


EC3.4.21- 


d6 


C-type lectin 7 


H. amigera 


gi|385202653 


34.99/5.98 


7.4 


ESTJHar 


128 


2 




d7 


C-type lectin 7 


H. amigera 


gi|385202653 


34.99/5.98 


5.5 


ESTJHar 


102 


2 




d8 


Serpin 


H. amigera 


gi|314912147 


43.53/5.04 


2.9 


ESTJHar 


165 


2 




dll 


Ferritin 


M. sexta 


gi|7272336 


26.51/5.44 


2.6 


ESTJHar 


193 


3 




Proteins decreased in diapause-destined larval hemolymph 


a2 


Growth-blocking peptide 
binding protein 


Heliosthis virescens 


gi|238915973 


23.63/5.16 


1.7 


EST_others 


79 


1 




a3 


Sulfa tase 


D. plexippus 


gi]3S7631 807 


56.65/5.34 


2.2 


EST_others 


105 


1 


EC3.1.6.14 


a4 


Sulfatase 


D. plexippus 


gi|357631807 


56.65/5.34 


2.3 


EST_others 


114 


1 


EC3.1.6.14 


a6 


EndoU protein, partial 


Popilio polytes 


gi|389615052 


28.70/8.38 


3.3 


ESTJHar 


217 


4 


EC3.1- 


a7 


Kazal-type inhibotor 


D. plexippus 


gi]357625979 


22.97/5.07 


1.8 


EST-Har 


229 


2 




b1 


Protein disulfide isomerase 


B. mori 


gi| 112984454 


55.59/4.60 


5.7 


EST_others 


140 


2 


EC5.3.4.1 


b3 


Ommochrome-binding 
protein 


M. sexta 


gi|400673 


30.61/5.89 


2.8 


ESTJHar 


207 


3 




b5 


Hemolymph proteinase 18 


M. sexta 


gi|56418417 


44.14/5.94 


2.0 


ESTJHar 


49 


1 


EC3.4.21- 


c2 


Serine proteinase-like 
protein 1 


H.armigera 


gi|208972549 


46.51/5.43 


1.8 


NCBI 


129 


3 


EC3.4.21- 


c3 


C-type lectin 


H. armigera 


gi| 106897087 


38.01/5.42 


1.9 


ESTJHar 


49 


1 




c4 


C-type lectin 


H. armigera 


gi| 106897087 


38.01/5.42 


2.9 


ESTJHar 


118 


1 




d1 


Serine proteinase-like 
protein 1 


H. armigera 


gi|208972549 


46.51/5.43 


10.7 


NCBI 


84 


2 


EC3.4.21- 


d2 


Serine proteinase-like 
protein 1 


H. armigera 


gi|208972549 


46.51/5.43 


10.7 


NCBI 


148 


3 


EC3.4.21- 


d12 


Unknown secreted protein 


Papilio xuthus 


gi|389609461 


1 7.50/8.55 


2.3 


ESTJHar 


195 


3 


EC3.4.21 



a Mowse scores greater than the significance threshold indicate identity or extensive homology (p < 0.05). Three database have their respective significance 
thresholds, the EST_Har database is 31, the NCBInr database is 46, and the NCBI EST_others database is 57. Additional data for single peptide-based identifications 
are presented in Additional file 2: Figure S2. 



at changed levels in the middle stage of the sixth instar: 9 
were found at higher concentrations, and 4 were detected 
at lower concentrations. y-Aminobutyric acid (GABA) 
was only detected in diapause-destined larvae, and other 8 
metabolites were present at higher concentrations in 
diapause-destined larvae: glycine, fructose, glucose, 
glycerol, ribitol, pyruvate, urea, and phosphoric acid. The 



concentrations of ornithine, 3-hydroxy-3-methylglutaric 
acid, phosphorylethanolamine, and ribonic acid were 
decreased. In the late stage of the sixth instar, all 4 
metabolites with altered profiles in diapause-destined 
larvae were found at higher concentrations: 7V-Acetyl- 
glucosamine, glucose, A/-Acetylglutamate, and gluconic 
acid. As shown in Figure 3A, more metabolites showed 
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Table 2 Metabolites identified in tha larval hemolymph 
from H.armigera 



Metabolite 


Fold change of relative ratio a 






Mid-late 


Early 


Mid 


Late 




5th 


6th 


6th 


6th 


Amino acid 


Alanine 


ns b 


ns 


ns 


ns 


p-Alanine 


ns 


ns 


ns 


ns 


Glutamate 


ns 


ns 


ns 


ns 


Glycine 


ns 


ns 


1.7 


ns 


Histidine 


ns 


-4.5 


c 


- 


Homoserine 


ns 


ns 


- 


- 


soleucine 


ns 


ns 


ns 


ns 


Leucine 


ns 


ns 


ns 


ns 


Methionine 


ns 


ns 


ns 


ns 


N-Acetylglutamate 


ns 


ns 




1.9 


Ornithine 


ns 


ns 


-2.4 


ns 


Phenylalanine 


ns 


ns 


ns 


ns 


Proline 


ns 


ns 


ns 


ns 


Serine 


ns 


ns 


ns 


ns 


Threonine 


ns 


ns 


ns 


ns 


Tryptophan 








ns 


Tyrosine 


ns 


ns 


ns 


ns 


Valine 


ns 


ns 


ns 


ns 


Sugar and polyols 


Arabinitol 


ns 


ns 






Cholesterol 


ns 


ns 


ns 


ns 


Erythritol 


ns 


ns 






Fructose 


ns 


ns 


3.8 


ns 


Galactitol 


ns 


ns 






Glucose 


ns 


1.7 


3.0 


3.1 


Glycerol 


ns 


ns 


2.2 


ns 


Myo-inositol 


ns 


ns 


ns 


ns 


Ribitol 


ns 


ns 


2.4 




Sitosterol 


ns 


ns 


ns 


ns 


Trehalose 


ns 


8.3 


ns 


ns 


Turanose 




ns 


ns 


ns 


Metabolic intermediates 


Fumaric acid 


ns 


ns 


ns 


ns 


Isocitric acid 


ns 


ns 


ns 


ns 


Malic acid 


ns 


ns 


ns 


ns 


Pyruvate 


ns 


ns 


1.5 


ns 


Succinic acid 


ns 


ns 


ns 


ns 


Urea 


ns 


ns 


2.5 




Small molecules 


Free amine 


ns 


ns 


ns 


ns 


Free phosphate 


ns 


ns 


ns 


ns 


Phosphoric acid 


ns 


ns 


1.5 


ns 



Table 2 Metabolites identified in tha larval hemolymph 
from H.armigera (Continued) 



Other metabolites 


2-Hydroxymandelic acid 


- 


ns 


ns 


3-Aminoisobutyric acid ns 


ns 


ns 


ns 


3-Hydroxy-3-methylglutaric 
acid 


- 


-2.7 


ns 


3-Hydroxypropanoic acid ns 


ns 


" 


" 


ns 


ns 


ns 


ns 


Galatonic acid ns 


ns 






Y-Aminobutyric acid ns 


ns 


D.specific 


ns 


Gluconic acid ns 


ns 


ns 


1.7 


Glyceric acid 


ns 






Linoleic acid 


ns 






Methylcitric acid 


ns 


ns 


ns 


Mevalonic acid ns 


1.7 




ns 


N-Acetylglucosamine ns 


ns 




5.3 


Octadecenoic 




ns 


ns 


Palmitic acid ns 


ns 


ns 


ns 


Pantothenic acid ns 




ns 


ns 


Phosphorylethanolamine 1.5 


ns 


-2.4 


ns 


Ribonic acid ns 


ns 


-2.9 


ns 


Stearic acid 




ns 


ns 


Threonic acid ns 


ns 







a Fold change with a positive number indicates metabolic levels were 
significantly increased in diapause-destine larval hemolymph, while a negative 
number indicates metabolite levels were significantly decreased in destined 
diapause-destine larval hemolymph. Based on four biological replicate analyses, 
only metabolites that changes >1 .5-fold in relative ratios (p-vaiue < 0.05} were 
considered for further analysis, and were marked in bold. b ns, no significant 
difference. c -, metabolite was found in less than three measurements and 
wasn't able to be quantified accurately. d Metabolite was only found in 
dispause-destined larval hemolymph. 

altered concentrations in the diapause preparation 
phase than in the diapause induction phase; this re- 
sult is consistent with the proteomic analysis above. 
Among the metabolites with altered profiles (Figure 3B), 
other metabolites accounted for 36% (8/22) of the total. 
Sugars and polyols represented 32% (7/22) of the total, 
and amino acids made up another 18% (4/22). 

The results of principal component analysis (PCA) 
carried out on the 59 identified metabolites suggest that 
there are no obvious physiological differences between 
the nondiapause- and diapause-destined fifth instar 
larvae (Figure 4A). However, physiological differences 
are observed in the sixth instar larvae (Figure 4B, 4C, 
and 4D) by PCI or PC2. 

Discussion 

Prediapause phase is a key period during which insects 
determine whether to develop into adults without inter- 
ruption or to enter diapause at the pupal stage. Taking the 
latter path requires changes in the individual's metabolic 
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Figure 3 Number (A) and classification (B) of significantly changed metabolites in H. armigera larval hemolymph in the 
prediapause phase. 



pathway to elevate energy reserves and produce cryopro- 
tectants to defend against low temperature. However, little 
is known about the relationship between environmental 
signals and the expression of the diapause phenotype [3]. 
A major reason for this lack of knowledge is the fact 
that no obvious phenotypic changes can be observed 
in diapause-destined individuals during the prediapause 
phase. In this study, we combined proteomic and metabo- 
lomic approaches to investigate changes in the expression 
of proteins and metabolites between larval hemolymph 
from nondiapause- and diapauses-destined H. armigera. 
Significantly changed proteins and metabolites were 
identified and analyzed, and our results revealed important 
clues regarding the diapause induction and preparation in 
the prediapause phase. It is well known that the short 
daylength is the major factor behind the effective induction 
of diapause, while low temperature is an associated minor 



factor. We believe that the altered levels of proteins and 
metabolites in hemolymph are the result of a cooperative 
influence of photoperiod and temperature because we used 
short daylength and low temperature in our experiment. 

Diapause induction 

Pyridoxine 5 -phosphate oxidase, an enzyme that is 
increased in the early stage of the sixth instar, catalyzes 
the biosynthesis of pyridoxal 5 -phosphate [11]. The vitamin 
B6 derivative pyridoxal 5 -phosphate acts as a coenzyme 
and participates in the biosynthesis of some neurotransmit- 
ters, such as serotonin, dopamine, and GABA [12]. Biogenic 
amine neurotransmitters in hemolymph, including 
serotonin and dopamine, have been shown to be involved 
in insect diapause in Mamestm brassicae and in Pieris 
brassicae [13,14]. In the middle stage of the sixth instar, 
GABA was only detected in diapause-destined larval 
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hemolymph. GABA is a major inhibitory neurotransmitter 
in insects. It has been reported that GABA and the GABA 
receptor antagonist picrotoxin influenced the incidence 
of diapause in Sarcophaga bullata [15]. GABA was 
also found to repress diapause hormone (DH) secretion, 
leading to diapause in Bombyx mori [16]. We speculate that 
the high concentration of GABA in diapauses-destined 
H, armigera larvae may be involved in inhibiting neurohor- 
mone secretion. This leads, in turn, to changes in larval 
metabolism and finally causes insects to enter diapause. 
Therefore, the increased pyridoxine 5 -phosphate oxidase 
in the diapause induction phase is most likely a response 
to the short daylength signal to regulate endogenous 
chemicals such as neurotransmitters for insect diapause. 

The amino acid histidine is a precursor to histamine, 
and the histidine levels were dramatically decreased in the 
early stage of the sixth instar. There is much evidence to 
suggest that histamine acts as a neurotransmitter released 
from the photoreceptor synapse of insects and other 
arthropods. Histamine acts by directly gating chloride 
channels and inhibits the activity of central nervous 
circuits in arthropods [17,18]. Diapause is mostly induced 



in H. armigera by short daylength signals which are 
received by photoreceptors in the brain. Significant 
changes in the levels of histidine most likely indicate 
altered levels of the photoreceptor transmitter histamine, 
which is likely to participate in the molecular events of the 
diapause induction phase. 

Metabolic storage 

Diapause-destined individuals must store additional meta- 
bolic reserves such as fuel reserves and cryoprotectant 
molecules because adequate reserves are critical to allowing 
individuals to sustain a basic metabolism in the harsh 
environmental conditions that occur during diapause. 
The major energy resource glucose was maintained at 
a constant high level during the sixth instar in diapause- 
destined larvae. Fructose and pyruvate levels were also 
elevated in the middle stage of the sixth instar. In our 
previous study, glucose and pyruvate levels decreased in 
diapausing pupae whose metabolic activities were dra- 
matically reduced [9]. Accordingly, accumulation of these 
glycolysis-related metabolites in diapause-destined larvae 
implies increased energy reserves in vivo. Interestingly, 
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trehalose, as an important cryoprotectant and energy 
resource that usually accumulates at the pupal stage [19], 
was heavily elevated in the early stage of the sixth instar. 

On the other hand, there is significant evidence to 
suggest that synthesis of cryoprotectants is associated with 
insect diapause [20]. Glycerol and ribitol are important 
cryoprotectant molecules found in insects [21]. These 
two cryoprotectants were found to accumulate in the 
diapausing flesh fly, Sarcophaga crassipalpis [22], and 
in the rice stem borer, Chilo suppressalis [23]. In this 
study, both compounds were found to be elevated in 
the middle stage of the sixth instar, showing that polyol 
concentrations in diapause individuals are closely 
correlated with cold resistance. Two major sugars, glucose 
and fructose, were also elevated in the sixth instar. 
These two sugars may also be used as precursors for 
the synthesis of cryoprotectants [24]. 

Urea levels were increased in diapause-destined larval 
hemolymph in the middle stage of the sixth instar. On 
the other hand, the levels of ornithine, which serves as 
precursor of urea, were decreased. The concentrations 
of 7V-Acetylglutamate, an activator of carbamoyl phosphate 
synthetase in the urea cycle [25], were also increased in 
diapause-destined larvae during the diapause preparation 
phase (the late sixth instar). It has been reported that the 
overwintering wood frog, Rana sylvatica, can accumulate 
additional urea, suggesting that the organic osmolyte 
also functions as a cryoprotectant such as glycerol in 
terrestrially hibernating amphibians [26]. Thus, urea in 
diapause-type insects may contribute to cold resistance. 

Although both nondiapause- and diapause-destined 
larvae were reared at the same temperature (20°C) with dif- 
ferent photoperiods, diapause-destined individuals could 
accumulate additional polyols and sugars. This result 
indicates that the synthesis of cryoprotectants and energy 
reserves is mainly attributable to short daylength conditions 
in this study, and the photoperiod is a key factor 
influencing the induction of diapause. Experiments in 
other insect species also showed this special linkage 
between photoperiod and cold-hardiness [20]. 

Immune response 

Insects are vulnerable to fungal and bacterial infections 
during diapause. The innate immune system continues 
to function in diapausing individuals and plays an 
important role in preventing the invasion of the body 
by microorganisms [27]. Serine proteases participate 
in many aspects of invertebrate immunity, including 
hemolymph coagulation, antimicrobial peptide synthesis, 
and melanin synthesis [28]. The activity of these enzymes 
is regulated by serine protease inhibitors (serpins). Thus, 
serpins in insect hemolymph also function to protect 
individuals from infection [29]. The expression of serine 
protease-like protein 2, serpin 3a, and serpin was found to 



be increased in diapause-destined individuals. On the 
other hand, levels of Kazal-type inhibitor, hemolymph 
proteinase 18 (a serine proteinase identified in Manduca 
sexta hemolymph), and serine proteinase-like protein 1 
were decreased in these individuals. Serine protease-like 
protein 2 and serpin, which accumulate at the late sixth 
instar of diapause-destined larvae, may be used in the long 
diapause phase to prevent invasion by microorganisms. 
This function has been reported for Nasonia vitripennis 
[30], Sesamia nonagrioides [7], Leptinotarsa decemlineata 
[31], Delia antiqua [32], and Culex pipiens [33], indicating 
that serine protease-like protein 2 and serpin are most 
likely closely linked to insect diapause. We speculate 
that the serine proteinase-like protein 1 may not be 
suitable to act in diapausing individuals at low tem- 
perature, and the expression of this protein is conse- 
quently decreased in the late sixth instar of diapause- 
destined larvae. 

Two different C-type lectins exhibited significantly altered 
expression in the diapause preparation phase. The 
expression of C-type lectin 7 was increased at the late stage 
of the sixth instar in diapause-destined larvae, while the 
expression of another C-type lectin was decreased. C-type 
lectins are calcium-dependent carbohydrate-binding 
pattern recognition proteins. Lepidopteran C-type lectins 
have been shown to function in innate immune responses 
such as phagocytosis, prophenoloxidase activation, and 
hemocyte nodule formation [34,35]. C-type lectin 7, which 
accumulates at the late stage of the sixth instar in diapause- 
destined larvae, may be highly stable at low temperature 
and can be used by diapausing individuals. Pattern 
recognition proteins can recognize foreign invading 
microorganisms, which is the crucial first step in producing 
an immune response [36]. A pattern recognition protein, 
beta-l,3-glucan recognition protein 3 ((3GRP3), was found 
to exhibit increased expression at the late stage of the 
sixth instar in diapause-destined larvae. Fungal infec- 
tion can be recognized by (3GRP through its binding 
to beta-l,3-glucan on the surface of fungal cell walls. 
Subsequently, (3GRP may work in concert with serine 
proteases to trigger the prophenoloxidase activation 
pathway, which is an important defense against microbial 
invasion in insects [37]. 

Previous studies on serine proteases, protease inhib- 
itors and lectins in hemolymph have revealed their 
functions in immune response for morphogenesis. Because 
both diapause- and nondiapause-type individuals must 
undergo a molting process, high expression levels of prote- 
ases, protease inhibitors, and lectins in these two types of 
individuals are correlated with morphogenesis. However, 
high expression levels of these proteins only in diapause- 
type individuals may serve in immune responses to 
facilitate survival in adverse environmental conditions 
over the long diapause period. 
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Stress tolerance 

One of the most common features in insect diapause 
is enhanced stress tolerance, which includes resistance 
to highly reactive chemical species, desiccation, ultraviolet 
radiation, salt, and other environmental or physio- 
logical stresses [38]. The expression of several stress 
tolerance-related proteins was significantly altered in 
the larval hemolymph of diapause-destined individuals. 
The expression of glutathione S-transferase (GST) was 
increased in the early stage of the sixth instar larvae. We 
showed that both GST mRNA and protein were heavily 
expressed in pupal brains of H. armigera at diapause 
initiation [6,39]. It has been reported that diapause 
hormone expressed at the larval stage of B. mori, as 
an intensified memory, can subsequently induce high 
expression of diapause hormone at the pupal stage 
[40]. Therefore, we speculate that GST expressed in 
the early stage of the sixth instar larvae may cause 
high GST expression at the pupal stage. GST has 
antioxidant activity and can protect cellular components 
of diapause individuals from oxidative damage, and thus, 
the accumulation of GST in diapause individuals is most 
likely involved in an oxidative stress response [41]. The 
levels of ferritin were also increased in the late stage of the 
sixth instar. High abundance of ferritin was also found in 
other diapause-destined insects such as N. vitripennis [30] 
and Calanus finmarchicus [42], and in the pupal brains of 
H. armigera at diapause initiation [39]. It has been shown 
that the concentration of ferritin increases in response 
to stresses such as anoxia [43]. Ferritin serves as an 
important ion storage protein in hemolymph and 
functions in maintaining the iron balance. The accumula- 
tion of ferritin can prevent destructive oxidative reactions 
by mitigating iron toxicity [44]. 

In addition, jV-Acetylglucosamine (GlcNAc) levels were 
found to be elevated in diapause-destined larvae in the late 
stage of the sixth instar. GlcNAc is the monomeric unit of 
the polymer chitin, which is one of the major components 
of insect cuticle [45]. Cuticle structural components are 
enriched in the diapause preparation phase to enhance 
waterproofing and to promote survival [46]. Furthermore, 
increased pupal cuticle protein was observed during early 
adult diapause in the mosquito G pipiens [47]. Elevated 
levels of GlcNAc most likely contribute to a harder cuticle 
in diapause individuals, improving tolerance to winter 
stresses and defenses against microbial invasion. 

Others 

Mevalonic acid was increased in the early stage of the sixth 
instar. Mevalonic acid is a precursor to the mevalonate 
pathway, which is the early step of juvenile hormone (JH) 
synthesis in insects [48]. JH participates in the regula- 
tion of insect development, metamorphosis, diapause, 
and reproduction [1,49,50]. To date, no convincing 



experimental evidence exists to support the notion 
that JH plays a role in the induction of pupal diapause. 
Although major pulses of JH activity with a rhythmicity of 
24 h were unique to pupal diapause-destined individuals 
in the flesh fly S. crassipalpis [51], the function of JH in 
the onset of pupal diapause requires further research. 

An imaginal disc growth factor-like protein showed 
significantly increased expression in the middle stage of 
the sixth instar. Imaginal disc growth factors (IDGFs) 
are soluble growth factors that promote cell proliferation 
in imaginal discs. Several IDGFs have been identified 
from the fat body of Drosophila; these growth factors are 
secreted and transported to target tissues through 
hemolymph [52]. Thus, the IDGFs that accumulate in 
diapause-destined larval hemolymph may be used for the 
development of imaginal discs when diapause is broken. 

A number of other proteins and metabolites showed 
significantly altered expression in diapause-destined 
larvae. These include sulfatase, endoU protein, protein 
disulfide isomerase, glycine, phosphorylethanolamine, 
and ribonic acid. Further research is required to reveal the 
roles of these proteins and metabolites in the regulation 
of diapause. 

Conclusions 

In summary, we identified significantly changed proteins 
and metabolites in larval hemolymph of diapause-destined 
H. armigera using parallel proteomic and metabolomic 
analyses. In the prediapause phase, external environmental 
signals are received and processed by larval brain, and 
subsequently produced endocrine signals induce changes 
of proteins and metabolites in hemolymph. Proteins and 
metabolites that altered in the prediapause phase were 
involved in various aspects of diapauses: diapause induc- 
tion, metabolic storage, immune response, stress tolerance, 
and others. These results here provided many valuable 
insights into the molecular characteristics of diapause 
induction and preparation. However, more research 
needs to be performed on the functions of altered 
proteins and metabolites in the future to understand 
the molecular mechanism in the prediapause phase. 

Methods 

Insects rearing and tissue collection 

H. armigera was kindly provided by Dr. J.-Y. Su (Nanjing 
Agriculture University, China). Larvae were reared on an 
artificial diet at 20°C and long daylength conditions 
(14 L : 10 D) to generate non-diapause pupae, and at 20°C 
and short daylength conditions (10 L : 14 D) to generate 
diapause pupae (Additional file 4: Figure S4). Diapause was 
determined by the retention of eyespots in the postgenal 
region of pupa [53]. Larval hemolymph was collected 
into microcentrifuge tubes containing phenylthiourea 
by clipping off a proleg with microscissors. Each sample 
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contained hemolymph pooled from 30 larvae. After 
centrifuged for 5 min at 1 2,000 g and 4°C, hemolymph 
samples were stored at -80°C until use. 

Protein extraction 

Three biological replicates were prepared for proteomic 
analyses. Protein extraction was first performed by 
precipitating the hemolymph samples in 10% TCA 
containing 20 mM dithiothreitol (DTT) on ice for 
1 h. After centrifuging at 1 2,000 g for 15 min, precipitate 
was washed three times with cold acetone (containing 
20 mM DTT). The pellet was dried and resolved in 
350 (iL rehydration buffer containing 7 M urea, 2 M 
thiourea, 4% (w/v) CHAPS, 40 mM DTT, 0.5% (v/v) 
IPG buffer (pH 4-7, linear), and 0.001% (w/v) bromophenol 
blue. Protein concentrations of the samples were 
determined by Bradford method [54]. 

Two-dimensional gel electrophoresis (2-DE) 

The hemolymph protein samples (340 \aL, containing 
approximately 1000 ^(g) were loaded onto 18 cm immo- 
biline pH gradients (IPG) strips (pH 4-7, linear; GE 
Healthcare), and applied by in-gel rehydration overnight 
in the reswelling tray. The first dimension isoelectric 
focusing (IEF) was carried out in an IPGphor II 
(GE Healthcare) as follows: 200 V for 1 h, 500 V for 1 h, 
1000 V for 2 h, gradient from 1000 V to 8000 V within 
1 h, and 8000 V for 5 h. After IEF, the IPG strips were 
incubated for 15 min in equilibration buffer containing 
6 M urea, 30% (w/v) glycerol, 2% (w/v) SDS, 50 mM 
Tris-HCl (pH 8.8), and 1% (w/v) DTT, and were incubated 
for another 15 min in the same equilibration buffer except 
that DTT was replaced with 2.5% (w/v) iodoacetamide. 
The equilibrated IPG strips were transferred to 12% 
SDS-polyacrylamide gels and sealed with agarose. The 
second dimensional electrophoresis was performed at 
80 V for 4 h, 100 V for 4 h, and 120 V for 2 h. The gels 
were stained with a mass spectrometry compatible 
silver staining method, without glutaraldehyde in the 
sensitization solution. Three independent biological 
replicates experiments were performed. 

2-DE image analysis 

The stained gels were scanned with an Image Scanner 
III (GE Healthcare). The gel images were analyzed with 
ImageMaster 2D Platinum 7.0 according to the user 
manual provided by GE Healthcare. Spots were detected 
and matched automatically; manual edition was carried 
out when necessary. For comparison of spot quantities, 
spot volumes were normalized as a relative volume 
(% vol, a percentage of the total volumes of all of the 
spots in the gel). Based on triplicate analyses, spots 
that showed a statistically significant (/?-value < 0.05, 



One-way ANOVA) and changed >1.5-fold were considered 
for further analysis. 

Mass spectrometry identification and database searching 

Spots were excised manually and subjected to in-gel 
digestion. Before digestion, spots were washed twice in 
Milli-Q water and dehydrated with acetonitrile, then 
reswelled on ice for about 45 min in 5 piL 25 mM 
NH4HCO3 digestion buffer containing 12.5 ng/uL 
porcine trypsin (Promega). Excess buffer that couldn't 
absorb by the gel spots was removed and replaced 
with digestion buffer without trypsin. Digestion was 
carried out at 37°C overnight. 

A sample (5 fiV) of each digest was applied directly to 
an AnchorChip plate (Bruker Daltonics), followed by 
1 piL matrix solution of 0.4 mg/mL oc-cyano-4-hydroxy- 
cinnamic acid (HCCA) in 70% acetonitrile, and then 
1 ftL TFA to remove remaining salt ions. Mass spectra 
were recorded on a Bruker Ultraflex III MALDI TOF/TOF 
mass spectrometer (Bruker Daltonics). Acquired MS and 
MS/MS spectra were sent as combined peak lists by 
BioTools software 3.1 (Bruker Daltonics) to Mascot search 
engine (Matrix Science), and searched against an in-house 
H. armigera EST database constructed in our laboratory 
(accession Nos. LIBEST_027843 NBR, LIBEST_027842 
DBR, LIBEST_027820 NFB and LIBEST 027821 DFB) [9], 
the NCBInr database or the EST others database. Only 
matches with significant scores (p < 0.05) were considered. 

The identified proteins were submitted to Universal 
Protein Knowledgebase (UniProt) to obtain KEGG 
annotations, and then they were classified according 
to Kyoto Encyclopedia of Genes and Genomes (KEGG) 
pathway maps. Some of them could not be assigned 
in KEGG pathway maps, therefore they were classi- 
fied according to the biological process of Gene 
Ontology. 

Metabolites extraction and derivatization 

Four biological replicates of samples were prepared. 
For each sample (contained 100 fA, hemolymph), 800 
methanol and 100 fiL Milli-Q water (containing 0.1 mg 
sucrose as an internal control) were added. The mixture 
was mixed thoroughly on a vortex and incubated on ice 
for 10 min. The mixture was vigorously extracted 
once again for 2 min and incubated on ice for 2 h. 
After centrifugation at 20,000 g for 10 min, 600 
supernatant was collected to dry in a vacuum concentrator. 
The dried sample was dissolved in 90 fiL methoxylamine 
hydrochloride (15 mg/ml in pyridine), incubated for 
16 h at room temperature, and trimethylsilylated in 
90 nL MSTFA/TCMS reagent (containing 1% TMCS) 
for 1 h. At last, 120 hexane was added to stop the 
derivatization reaction. 
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GC-MS analysis 

One microliter of the derivatized sample was auto-injected 
at 1:10 split ratio into a Trace GC Ultra-DSQII GC-MS 
(Thermo) equipped with a DB-5msUI column (length 
30 m, I.D. 25 mm, Agilent). Mass spectra were recorded 
from 50 to 450 m/z. The injection temperature was 280°C. 
The oven temperature was increased from 50 to 300°C at 
5°C/min, and maintained at 300°C for 5 min. 

Metabolites were at first identified by comparison of 
their retention time with authentic standards. Further 
identifications were achieved through spectral matching 
against the National Institute of Standard and Technology 
(NIST). Matches with a reverse similarity index (RSI) 
value >700 were considered. But some of the matches 
could not be linked to a clearly defined biological function, 
and therefore these metabolites were considered unidenti- 
fied. Individual integrated peak areas of the metabolites 
were all converted to response ratios relative to the internal 
standard sucrose. According to four biological replicates, 
metabolites that changed > 1.5-fold in relative ratios with a 
j?-value < 0.05 (One-way ANOVA) were considered for 
further analysis. Principal component analysis (PCA) 
was carried out with MetaGeneAlyse v 1.7.1 (Max Planck 
Institute of Molecular Plant Physiology, Germany). 
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